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Abstract 

We present a scaled particle density functional study of two-dimensional binary mixtures of hard 
convex particles with one or both species being ellipses. In particular, we divide our study into 
two parts. The first part is devoted to the calculation of phase diagrams of mixtures with the 
same elliptical shapes, but with (i) different aspect ratios and equal particle areas, (ii) equal aspect 
ratios and different particle areas and (iii) with the former and the later being different. In the 
second study we obtain the phase diagrams corresponding to crossed-mixtures of particles with 
species having elliptical and rectangular shapes. We compare the phase diagram topologies and 
explain their main features from the entropic nature of particle interactions directly related to the 
anisotropies, areas, and shapes of species. The results obtained can be corroborated by experiments 
on granular rods where the entropic forces are very important in the stabilization of liquid-crystal 
textures at the stationary states. 
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I. INTRODUCTION 



The demixing transition in mixtures of hard convex three-dimensional particles have 
received a considerably attention from the theoretical point of view [l- 14]. This is due to 
the possibility of formulating relatively simple theoretical models for mixtures of anisotropic 
particles which are not allowed to overlap. These models are based on the simple Onsager 



second- virial approximation 



15] or its rescaled version, the Parsons approach 16(. Also, the 



scaled particle theory (SPTV proposed firstly for hard spheres [17j] and later extended to 
other particles geometries |18l-l22|. has been a fruitful theoretical tool in the study of these 
kind of mixtures. All these models predict that a binary mixture of asymmetric enough 
particles can demix into two different phases, each one rich in different species. The nature 
of uniform demixed phases can be isotropic (I), a fluid with randomly oriented particle 
axis, or nematic (N) in which the particles' axes are oriented in average along the nematic 
director. All possible demixing transitions, i.e. I-I, N-N and I-N were found from the 
numerical implementation of these models. 

In two dimensions the theoretical works devoted to the study of the demixing transitions 
are more scarce, because strictly two-dimensional fluids occur less often in nature compared 
to the three-dimensional fluids. Examples of the former constitute the quasi-two-dimensional 
films obtained from preferential chemical adsorption of certain species at interfaces. The ad- 
sorption of polystyrene spheres [23|] or ellipsoidal latex particles |24J into a two-dimensional 
interface has allowed the experimental study of freezing and percolation/jamming transitions 
in two dimensions respectively. Also, the time-dependent behavior of anisotropic colloids 
near interfaces was recently studied through a theoretical modelling of the Langevin dy- 
namics of adsorbed particles at a flat interface |25lj . For a recent review on the behavior of 
colloidal particles at liquid interfaces see Ref. [261 ] . 

However, the most paradigmatic example constitutes the collection of granular anisotropic 
particles confined between two parallel plates with the distance between them less than 
the particle widths. A growing number of experimental works conducted on granular fluids 
subjected to different vibrational modes show spatial pattern formations, demixing scenarios 
of granular species and collective motions j^]. Those systems in which granular rods are 
vertically shacked have shown that the entropic nature of interactions between particles is 
very important in the explanation of phenomenology found. In particular, the presence of 
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defects of the nematic director due to the difference between particle orientations near and 



far from the walls was found 



28j |. Recent experiments on granular rods reported the presence 



of certain liquid- crystal textures with N, tetratic (T), and smectic (Sm) symmetries as the 
stationary states of these systems j^]. The presence of strong dissipation due to the inelastic 
interactions between grains is an important ingredient to take into account in any theoretical 
model of a shaken granular systems. However, the presence of demixed stationary states of 
a granular mixture could be described by models based of density functional theory. For 
example the ranges of molar and packing fractions for which the demixing is observed could 
be approximately given by these models. 

The SPT shows that a two-dimensional mixture of hard convex bodies never demixes 
However, recently we have shown that the same theory can predict 



into two I phases 



30]. 



I-N and N-N demixing 
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321 ] . The only difference with respect to three dimensions is the 



second order character of the I-N transition at pressures bellow the tricritical point [3ll . [32] . 
The main purpose of the present study is to implement the SPT for (i) a mixture of hard 
ellipses (HEs) with different aspect ratios and/or particle areas and (ii) for crossed- mixtures 
composed of HEs and hard rectangles (HR). The most convex particle with a fixed aspect 
ratio that can be imagined in two dimensions is just an ellipse. This property in turn might 
be reflected in the phase behavior of mixtures composed of HEs. We will show that this 
is indeed the case by studying the phase diagram topology as a function of particle aspect 
ratio and particle area. The one-component fluid of HEs has been studied via Monte Carlo 
simulations. This fluid exhibits, for high enough aspect ratios, an I-N transition occurring 
via the disclination unbinding mechanism, while for higher packing fractions a first order 



fluid-solid transition takes place 33J]. It has been shown that the equation of state for the 



one-component fluid of HEs following the SPT compares reasonable well with the simulation 



results 



34]. 



To implement the present study we firstly find an analytical expression for the excluded 
area between two HEs with different axes {a^, bj} (i = 1,2). Zheng and Palffy-Muhoray 3^] 
developed a method based on the calculation of the distance of closest approach between 
an ellipse and a circle. Thus, after rescaling the circle they obtained the distance of closest 
approach between two ellipses. Finally, the excluded area can be obtained by the angular 
integration of the distance of closest approach. However, this integral should be evaluated 
numerically. We will follow another procedure, in which the excluded area between two dif- 
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ferent ellipses can be calculated from the parametrization of the distance of closest approach 
as a function of the angle that forms the point of contact between two ellipses with a fixed 
reference axis. The resulting excluded area thus depends on the complete elliptic integral 
with its parameter being a function of the ellipse eccentricities and on the relative angle 
between particles. Further, we also calculate the excluded area between an ellipse and a 
rectangle. The details of both calculations are relegated to the appendices. 

This paper is organized as follows. In Sec. HI] we present the theoretical model used for 
the present study: the SPT. In Sec. IHII we analytically derive the spinodal curves of the I-N 
second-order transition, i.e. the packing fraction as a function of composition for general 
mixtures of two-dimensional hard particles . Sec. II VI is devoted to the analytical derivation 
of the N-N spinodal curve with the constraint of parallel particle alignment. In Sec. M we 
present the phase diagrams of mixtures of HEs (Sec. IV Al) and HEs and HRs (Sec. IV Bp 
obtained from the numerical minimization of the free-energy functional. Finally, in Sec. 
IVII the main conclusions are drawn. The appendices contain the details of the derivation 
of the excluded area between two HEs (Sec. and HEs and HRs (Sec. [Bj and different 
mathematical aspects of the functional minimization procedure (Sec. [C]). 



II. SCALED PARTICLE THEORY 



The central quantity of any theoretical model which accounts for two-dimensional hard 
particle interactions is the pair excluded area to particle v due to the presence of particle \x 
with the constraint of a fixed angle <ft between their main axis. Because the hard core nature 
of interactions, both particles can not overlap and thus the excluded area A^ip) defines a 
body with a geometry directly related with those of particles. The present paper is devoted 
to the study of phase diagrams of mixtures composed of HEs with different characteristic 
lengths [different semi-axis (a M , b M )] and a crossed-mixture of HEs and HRs. In Appendices 
IA1 and |B] we obtain the analytical formulas for the excluded areas of these mixtures. Once 
the excluded area is given we can construct a scaled particle theory (SPT) for a general 



mixture 



17] . The excess free energy per particle following the SPT can be calculated to be 



^ = - in(i - V ) + ^^^((4°))), (1) 
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where p is the total number density and i] = p £\ a^aj is the total packing fraction with 
Xi the molar fraction of species % and <Zj their particle areas. We have defined the function 
A$(<f>) = (A^cfj) — a M — a u ) /2, and used the short-hand notation 

((4°j))= r dfahM r dfahvfflAMfa-fa), (2) 

'o Jo 



to define the double angular average of A^} (</>) with respect to the orientational distribution 
functions h^tfi). Note that the second order density expansion of Equation §1§ give us the 
usual second virial approximation 



(b 2 ) _ P 



2 



-J2 X ^ X A( A ^}}- (3) 



The fact that the coefficients ((A$)), instead of the excluded areas, enter into the SPT 
free-energy is due to the form of the derived SPT: through a Taylor expansion of the excess 
chemical potential corresponding to the scaled particle, evaluated at low densities. The 
upper limit of integration in Equation (J2J) is 7r due to the head-tail symmetry of particles, 
i.e. = h^Tr — <p). These functions measure the degree of orientation of the particles 

along the fixed nematic director. In the present paper we have selected a variational family 
for these functions as 

exp(A M cos 20) 

Kw) = , 7T \ , (4) 

with I n (x) the n-th order modified Bessel function. The perfect nematic alignment is reached 
when A M — > oo, which gives us h^((p) = 5(<p) (the Dirac-delta function), while the isotropic 
fluid is obtained in the limit A M — > for which h^(<f)) = tt^ 1 . Note that this parametrisation 
excludes the possibility of study of the tetratic phase with fourfold angular symmetry h fJi ((f)+ 
tt/2) = h^<f>). 

The ideal part of the free-energy per particle can be computed as 

<p id = In T] - 1 + ^ x„ jln x n + J #fyu(^) In [vr/i M (0)] j , (5) 

= lnr ? -l + ^x / ,|lnx At + A M ^|-ln/ (A / ,)|. (6) 

where we have inserted (jl]) into fl5]) to obtain ([6]). 

The Gibbs free energy per particle g = if + P '/ p (with ip = ip id + ip exc the total free-energy 
per particle) calculated at a constant fluid pressure P should be minimized with respect 



to the parameters A M (/i = 1,2) to obtain their equilibrium values. The expression for the 
pressure following the SPT is 

In Appendix O we write explicit expressions for the equations to be solved numerically to 
obtain the equilibrium values of A M . The common tangent construction of the function g(x), 
with x — Xi, allows us to calculate the coexisting values for x (and thus of r\ once the fluid 
pressure is fixed) at the demixing transition. 



III. I-N BIFURCATION ANALYSIS 



In this section we obtain an analytic expression for the I-N spinodal of a general mixture. 
When A,, C 1 we can approximate the orientational distribution functions h^cft) up to first 
order in A M by 

1 
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[1 + Xf, cos 20) . 



(8) 



Inserting this expression in the free-energy per particle p = pid + p exc we obtain, up to 
second order in A„ 



fj, fa/ 

with y = p/(l — T]) and where we have defined the coefficients 

1 



a 
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d<f)cos(2i<f))A$((f)), i = 0,l. 



The minimization of fl9]) with respect to A M give us 

dp _ Xf, 



OX, 



0. 



(9) 



(10) 



(11) 



^ju/^vXi, 

V 

The term of Eq. fill I) enclosed by square brackets can be put in the following matrix form 

HX = (I + 2yU) X = 0, (12) 

where A = (Ai, A2) T 7^ 0, I is the 2x2 identity matrix while U is the matrix with elements 
oc$x v . Equation (TT2T) has a nontrivial solution only if det (H) = 0, which gives us the 
following result: 



1 + 2yr + {2yf Xl x 2 I o$a$ - 



0. 



(13) 



where we have defined the coefficient 



r 



) 



(14) 



Taking into account that the expression enclosed by brackets in Equation (1T3]) is equal to 
zero for any mixture of convex bodies, we find that the packing fraction at the I-N transition 
can be calculated as 

rf = r(a > y* = -J-, (15) 
' l + j/*(a)' 2r' v ; 

where we have used the shorthand notation (a) = ^ x^a^. Following the SPT the pressure 

of the isotropic fluid at the bifurcation point can be calculated as 















- A 4 





(16) 



where is the perimeter of species p. 



IV. N-N DEMIXING 



In this section we obtain analytically the criterion for the existence of N-N demixing 
and its spinodal curve with the constraint of parallel particle alignment. This constraint is 
necessary to obtain an explicit analytic function for this spinodal, while the freely rotating 
case can only be solved numerically. The spinodal instability of the binary mixture with 
respect to phase separation can be computed by requiring that T-L{r],x) = where we have 
defined 



H{r],x) 



<9 2 $ 



1 +P2 



P\P2 



d 2 $ 

dpidp 2 



(17) 



with pi = Xip, while $ = py3 ex is the excess part of the free energy density. From Equations 
(0Q) and ffTTj) we obtain 

Hiv, x) = (1 + y ) 2 {1 + 2y - 4x(l - x)y 2 U} , (18) 

where we have defined y = 77/(1 — 77), and U = (q 2 — 1) aia 2 / (a) 2 , with q = A± 2 / v /a i a 2- 
Solving "H(?7,x) = for r/, we find the explicit expression 

1 



V 
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(19) 



The packing fraction and composition at the critical point are i]* = 1/q, and x* = a 2 /{ai + 
a 2 ), respectively. Finally, the pressure at this point is 



We conclude that demixing transition occurs only if U > 0, i.e. for q > 1. For mixtures of 
HEs we find q ee = l^fsE [1 — s~ 2 \ /ir > 1, where we have defined the coefficient 



with Hi the aspect ratio of particle i while E[x] is the complete elliptic integral of the second 
kind. The same coefficients for a mixture of HRs and for a crossed-mixture of HEs and 
HRs are q„ = (y/s + l/\fs) /2 > 1 and q eT = (y/s + 1/a/s) I \pH > 2/a/tt > 1 respectively. 
We see that for s — 1 (no matter what the ratio a 2 /ai is) no demixing occurs in mixtures 
of particles with the same shape, while the crossed-mixture demixes even in this case. In 
Figure [T] we plot the coefficients q a $ as a function of s. It is interesting to note that at 
s ~ 5.759 two of the curves intersect (those corresponding to the mixture of HEs and the 
crossed-mixture). To illustrate the effect of this cross-over on the phase behavior we plot in 
Figure [2] the spinodals of both types of mixtures for s = 2.5 and s — 8, which are located at 
different sides of the cross-over point. 

V. RESULTS 

All the phase diagrams we present in this section were calculated via the bifurcation 
analysis (see Section IIIII) for the case of continuous phase transitions (as the second order 
I-N transition) and from the numerical minimisation of the parametrised Gibbs free energy 
per particle (see Section [TT] and Appendix E]) for the case of I-N or N-N demixing transitions. 
The later procedure allows us to obtain the equilibrium orientational distribution functions 
hfj,(4>) at coexistence and the coexisting values of the composition and packing fraction of 
the demixed phases at a given pressure. We begin the presentation of the results from the 
most simple binary mixture: a mixture of HEs. Note that the hard disk (HD) is a special 
case of ellipse with equal semi-axis. 




(20) 



S = K 2 /Ki > 1, 



(21) 



S 
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FIG. 1: The coefficient q as a function of s = K2/K1 for a mixture of HEs (solid curve), for the 
crossed- mixture of HEs and HRs (dashed curve), and for a mixture of HRs (dotted curve). The 
open circle shows the cross-over point. 
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FIG. 2: Demixing spinodals of mixtures of HEs (solid line) and of the crossed-mixture of HEs and 
HRs (dashed line) for s = 2.5 (a) and s = 8 (b). 
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FIG. 3: Phase diagram of the binary mixture of HEs (species 1) with k± = 10 and HDs (species 
2). Both species have the same particle areas a\ = a%. The regions of I and N stability are 
correspondingly labeled. The I-N spinodal is plotted with dashed line while the demixing binodals 
are shown with solid lines. The grey region constitutes the zone of mixture instability. The 
tricritical point is shown with a filled circle. 

A. Mixtures of HE 

The phase diagram topologies we show here are similar to those obtained for mixtures of 



HRs and hard discorectangles (HDRs) in recent studies [3l|, |22j. These studies have shown 
all possible demixing scenarios that mixtures of two-dimensional hard convex particles can 
exhibit. However, we will show that the ellipse, being the geometry with the most convex 
anisotropy we can imagine, considerably enhance the demixing transition in mixtures of 
particles with elliptical shapes. 

The first binary mixture we study is composed of HEs and HDs. The former has an 
aspect ratio of K\ = 10 while the area of both species are the same (ai = a%). The resulting 
phase diagram is shown in Figure [3] as a pressure-HE molar fraction (xi) plot. As we can see 
from this figure, at low pressures the mixture exhibits a second order I-N phase transition up 
to pressure value of 0Pai ~ 3.7 (the location of a tricritical point) above which the system 
demixes into an I phase rich in HDs and a N phase rich in HEs. The demixing gap is wider as 
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FIG. 4: Phase diagram of two binary mixtures of HEs with K2) = (6, 2) (a) and («i, K2) = (5, 2) 
(b) while the particles areas of both species are equals. The meanings of different lines and labels 
coincide with those of Fig. [3l The tricritical and critical endpoints are shown with filled circle and 
square, respectively. Finally, the critical point is shown with an open circle. 

the pressure increases from the tricritical point. We should point out that the second order 
I-N transition obtained from density functional calculations replace the Kosterlitz-Thouless 
I-N transition obtained by simulations and theoretical models that properly account for the 
nematic director fluctuations. The SPT theory does not includes these fluctuations, thus 
predicting a continuous increase of the long range nematic ordering. 

The second study we have carried out concerns the calculation of the phase diagram of 
two similar binary mixtures of HEs, again with the same particle area ai = 02 but with 
different aspect ratios given by K\ — 6 and K2 = 2 for the first mixture [Figure H] (a)] while 
they are K\ — 5 and K2 = 2 for the second mixture [Figure |4] (b)]. There is an important 
difference between both phase diagrams in that the first one has a region of I-N demixing 
above the tricritical point [see Figure S] (a)] which is replaced at higher pressures by a N-N 
demixing (just when the I-N spinodal intersects the I binodal of the I-N coexistence). We call 
this point a critical end point. The I phase as usual is rich in HE species with small aspect 
ratio. The second phase diagram exhibits a second order I-N transition for all compositions 
and at higher pressures the system demix into two different N phases. This N-N demixing 
transition ends at a critical point. 
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X l 

FIG. 5: Phase diagram of a binary mixture of HEs with aspect ratios K\ = 5 and k<i = 2 while the 
partcle area a\ = 2.5a2- The meanings of labels and lines are as in Figigure[3j 

How affects the difference in particle areas on the phase behavior constitutes the following 
point we have elucidated. The aspect ratios of both ellipses are again k% = 5 and k-i = 2 
while the ratio between particle areas is a\/a2 = 2.5. In Fig. |5]we show the resulting phase 
diagram. Again appears a window of I-N demixing but now the stability region of the second 
nematic (N 2 ) shrinks considerably. 

Now we consider a mixture of HEs with the same aspect ratios K\ = k 2 = 5 and different 
areas, specifically a\/a,2 = 50. The phase diagram is shown in FigureEJ This mixture exhibits 
two tricritical points defining the lower and upper limits of the I-N demixing and two critical 
points defining the proper limits of the N-N phase separation. In Figure [6] the molar fraction 
of the first species is substituted by its area fraction, i.e. x a = X\a\j \x\a\ + X2&2)- I n such 
a way the position of all these points can be properly discerned. 

B. Mixtures of HE and HR 

This section is devoted to the study of mixtures of particles with different shapes, specif- 
ically HE- and HR-shaped particles. The first crossed-mixture we have studied is composed 
of species with aspect ratios coinciding with those of the already studied HE mixture with 
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FIG. 6: Mixture of HEs with the same aspect ratios K\ = ft 2 = 5 and different particle areas 
a\ = 50a2. The critical and tricritical points are shown with open and filled circles respectively. 

phase diagram plotted in Figure H] (b). Specifically the HR has the larger aspect ratio 
(«i = 5) while the HE has k 2 = 2, and the areas of booth particles coincide. The aspect 
ratio of the rectangular particles is selected so as to be large enough to forbid the T phase 
stability 3JJ. In Figure [7] we show the phase diagram of this mixture. In the same figure we 
plot for a better comparison the phase diagram of Figure H] (b) . We can extract two main 
differences between both phase diagrams. Firstly, the I-N transition occurs at lower pres- 
sures for the HE mixture. This feature is related to the change in the excluded the 
particle aligns preferentially along the nematic director. This change is larger in mixtures 
of HEs compared with crossed-mixtures. The other main difference is related to the size of 
the N-N demixing gap. The crossed-mixture has the wider gap which can be understood in 
terms of packing of particles at high pressures: particles with the same shape pack better 
than those with different shapes. This in turn has already been quantified in Section IIVI 
through the coefficients q ee and q er (see Figured]). 

The last study we have carried out is related to the effect that different particle shapes 
have on the phase diagrams of binary mixtures with species having very dissimilar areas 
(specifically a\ = 50ci2) but equal aspect ratios (which were fixed to 5). To this purpose we 
have calculated the phase diagrams of a mixture of HEs (the dotted lines in Figure [8]), a 
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FIG. 7: A mixture of HRs (species 1) and HEs (species 2) with the same particle areas a\ = (12 while 
the aspect ratios are K\ = 5 and k<i = 2, respectively. The dashed line represents the I-N spinodal 
while the solid lines represent the N1-N2 demixing transition. Also plotted is the phase diagram of 
Figure 0] (b) corresponding to the mixture of HEs with the same particle aspect ratios and areas. 
The dotted and dashed-dotted lines represent the spinodal and demixing curves, respectively, for 
this case. 

crossed-mixture of HEis and HRs with rectangles being the large species (the solid line in 
Figure [S} and the opposite mixture with the ellipses being the large species (the dot-dashed 
line in Figure [S]). As we can see from this figure the lower I-N spinodal corresponds to 
the mixture of HEs which is related to the major gain in the excluded area by the particle 
alignment of the most convex species, the HEs. For the same reason, we expect that the I-N 
demixing region is widened for crossed-mixtures with ellipses being the large species. The 
mixture of HEs and HRs, the later being the large species, has the smallest demixing gap 
which is related to the fact that a great amount of large rectangles enhances the stability 
of the isotropic phase with respect to the nematic, due to the disordered effect that cause 
the rectangular corners. Finally, the scenario with an intermediate demixing gap occurs for 
mixtures of HEs. 
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FIG. 8: Phase diagrams of the binary mixtures with equal particle aspect ratios (k\ = K2 = 5) but 
different areas (ai = 50a2). The dotted lines represent the mixture of HEs, the solid lines represent 
the mixture of HRs (species 1) and HEs (species 2) and the dot-dashed lines represent the mixture 
of HEs (species 1) and HRs (species 2). The positions of critical points are shown with different 
symbols. 

VI. CONCLUSIONS 

One of the main purposes of the present work has been to encourage the statistical me- 
chanical study of two-dimensional mixtures of hard anisotropic particles as a starting point 
to understand the phenomenology inherent in mixtures of granular particles. While three- 
dimensional mixtures have been systematically studied via density functional calculations 
and by MC simulations, there is a lack of results on two-dimensional mixtures, with scarce 
theoretical works and no simulation results. 

We have shown here that the phenomenology inherent to these mixtures are as rich as the 
three-dimensional mixtures, with all possible demixing scenarios where the demixed phases 
have isotropic or nematic symmetries. The absence of I-I demixing and the presence of second 
order I-N transition constitute the main differences with respect to the three-dimensional 
case. 

While mixtures of hard particles with rectangular and disco-rectangular geometries have 
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been recently studied with density functional tools (3l|, |32|], the present work constitutes the 
first attempt to incorporate the elliptical geometry into these studies. For a given aspect 
ratio the ellipse is the most convex particle we can imagine. As we have shown here, this 
geometrical property is translated to the the phase diagram topologies of mixtures with 
one or both species being ellipses. We trust that the same kind of results can be found in 
mixtures of ellipsoidal and cylindrical granular rods, the Basmati rice being an important 
paradigm of the former. 

We should take the conclusions about the stability of uniform phases that the SPT 
predicts with certain caution. Phase transitions to non-uniform phases, such as smectic or 
crystalline phases, can occur at pressures bellow the maximum values depicted in the phase 
diagrams plotted here. Only a full minimization of a density functional with respect to the 
density profile p(r, <ft) (depending also on the spatial variables) can clarify this problem. A 
density functional for mixtures of HEs based on the Parsons approach can in principle be 
implemented for the general density profile inhomogeneities. However this functional treats 
the translational degrees of freedoms in a crude approximation, wich is a serious drawback 
for two-dimensional systems where the many-body correlations are very important. 
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Appendix A: Calculation of the excluded area between two ellipses 

In this section we calculate the excluded area between two different ellipses. The center 
of an ellipse with semi-axis ai and b>i will be placed at the origin of the coordinate system 
with the major semi-axis bi being parallel to the y Cartesian axis. The other ellipse, with 
semi-axes a2, b2, is placed at the point ro = (xq, yo) and has its major axis forming and angle 
(f) with respect to the long axis of the first ellipse. The vector ro is such that there is only 
one common tangent point between both ellipses (see Figure [9]). 
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FIG. 9: Configuration of the closest approach between two ellipses with semi-axes {ai,fc>i} and 
{a2, ^2} with a relative orientation defined by the angle between their long axis, which is equal to 
0- 

The parametrised coordinates of the first ellipse are 

ri(*i) = M*i),?/i(fi)] = (ai cos *i,bi sin ^e[0,27r], (AI) 

while the second vector r 2 (t 2 ) = [^2(^2), 2/2(^2)] is parametrically defined as 

^2(^2) = xo + a 2 cost 2 cos0 — b 2 sint2sin0, (A2) 
2/2(^2) — 2/0 + b 2 sint 2 cos0 + a 2 cost 2 sin0 (A3) 

with t 2 G [0,27r]. The tangent unit vectors Tj corresponding to both ellipses, calculated by 
differentiating [xi(ti),yi(ti)] with respect to the parameters tj, are 

— — {— a i sinti, bi costi), (A4) 

n 

^"2(^2) = — (—32 sin t 2 cos — b2 cos £2 sin 0, b2 cos £2 cos — a 2 sin t 2 sin 0), (A5) 

?~2 

where Tj = \/ af sin 2 tj + bf cos 2 tj. The condition of a common tangent point between both 
ellipses can be expressed by the following set of equations: 

ri(ti) = r 2 (t 2 ), (A6) 
Ti(ti) = -r 2 (t 2 ), (A7) 
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which can be solved to obtain sint 2 (or cos £2) a s a function of t% and also for Xo(ti) and 
Uo(ti) with the result 

b2(bi costi sin0 — ai sinti cos0) . , 

smt 2 = , (A8) 

— a 9 (ai sintx sin0 + bi costi cos0) . . 

cost 2 = j^j- 2 , (A9) 

%o(ti) = ai costi — a 2 cost2 cos0 + b 2 sint 2 sin 0, (A10) 

2/o = bi sinti — b 2 sint 2 cos 4> — a 2 cost 2 sin0, (All) 

where 

5(ti) = a 2 (a 2 sin 2 1\ + b\ cos 2 1\) + (b 2 — a 2 )(bi costi sin — ai sinti cos 0) 2 

= b 2 (a 2 sin 2 ti + b 2 cos 2 ti) - (b 2 - a 2 )^ sinti sin0 + bi costi cos0) 2 , (A12) 

The excluded area A(4>) can be calculated as 

A ^) = \f [y'oitdMtd-yoitiKit^dh, (A13) 

After some long calculations we arrive at the following relatively simple formula for this 
function: 

A((f>) = Trfcib! + a 2 b 2 ) + 2 [A + (0) + A_(0)] (A14) 

where -E(k) = / \/l — « sin 2 0cZ0 is the complete elliptic integral of the second kind and 
where we have defined 



4A + (0)A_(0) 
[A + (0) + A_(0)f 



= r rs™r7L > (A15) 



A±(0) = v /(b 1 b 2 ±a 1 a 2 ) 2 -(b 2 -a 2 )(b 2 -a 2 )cos 2 0- (A16) 
Appendix B: Ellipse-rectangle excluded area 

The purpose of this section is to obtain an expression for the excluded area between an 
ellipse with semi-axes a and b (b > a) and a rectangle with length and width equal to L and 
a, respectively. We show in Figure H] an schematic representation of the excluded area. It 
can be seen from the figure that the total area can be computed by summing two different 
contributions: (i) one coming from four triangles (two of them having different areas Af\ 
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FIG. 10: The excluded area between an ellipse with semi-axes a and b and a rectangle of width a 
and length L. Different sections of the excluded area are shown. 



i = 1,2) and (ii) the second, coming from four figures with elliptical borders (two of them 
having Af\ i = 1,2). The corresponding to the triangles are 



where are the absolute values of the vectors shown in Figure HI These areas can be 
computed as 



To compute the areas of figures with the elliptical borders we define the radius vector u 
from the origin to the border of the excluded area (see Figure [TU]) . which in turns is equal to 
the sum of two vectors: the vector c = (c x , c y ) which accounts for the position of the center 
of mass of the rectangle with its inner vertex touching the perimeter of the ellipse and the 
vector r = (a cos#, bsin^) which has been parametrised with the angle 9 G [#i,6> 2 ]. Thus, 
the area of this figure can be calculated from 




(Bl) 



(B2) 





(B4) 




(B5) 
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with ip the angle that the vector u forms with the x-axis. Taking into account the fact that 
u = a/ (a cos 9 + c x ) 2 + (b sin 9 + c y ) 2 and 



sin ip = u 1 [c y + b sin 9) , cos ip = u 1 (c x + a cos ( 



(B6) 



we find from (1B5|) that 



A 



ab ( Of - 9f )+bc x [ sin 6%> - sin Bf 1 ) - ac^ ( cos 9 [ ?> - cos 0}' 



3 (i) 



"1/ 



3 (i) 



i(0 



(B7) 



where are the components of the vector describing the movement of the center of 
mass of the rectangle inside the two different regions with elliptical borders, while o\ % are 
the lower and upper limits of the parameter 9 inside these regions. All these quantities are 
explicit functions of a, b, a and L. It can be easily shown that 



(e) 



7T 



ab H — \/a 2 cos 2 6 + b 2 sin 2 d> H — \/a 2 sin 2 d> + b 2 cos 2 



(B8) 



Thus, the total excluded area can be obtained by summing all the contributions, i.e. A((f>) 
2£, 



A w + A (e) 



with the final result 



A{4>) = La + vrab + 2 



a\ a 2 sin 2 6 + b 2 cos 2 + L\ a 2 cos 2 + b 2 sin 2 



(B9) 



Appendix C: Minimization of the Gibbs free-energy per particle 

As we have already pointed out before, the double angular average is the mam 

ingredient of SPT which can be computed as 



(« ) ))= / WM^M), 



(Cl) 



where the functions 



= / d<f ) 'h II (<p')h u (<p + <f ) '), 



(C2) 



were defined, which using Equation (j3j) become 



7rJ (A At )/ (A I /) : 



£M) = x /A 2 +A 2 + 2A M A,cos2(/. 



(C3) 
(C4) 
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To numerically implement the minimization of the Gibbs free energy per particle g with 
respect to the parameters A M we have used the following expressions for the derivative of g 
with respect to A M 



dg 
dX„ 



x,, < A,, 



/o 2 (A.) 



MM 



y 



jJ.v ( -i 



dX„ 



(C5) 
(C6) 



where 



d%M = 1 
OX, 7r/ (A,)/ 2 (A M ) 



(A M + X u cos 2(f)) 



hit 



Hi/ j 



/o(A M )-/o(^)A(A M ) 



(C7) 



We have solved the set of two non-linear equations dg/dX^ = (/x = 1,2) with respect to 
two unknowns (Ai,A2). 
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